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Abstract. We consider two natural problems arising in geometry which are equivalent to 
the local solvability of specific equations of Monge- Ampere type. These are: the problem of 
locally prescribed Gaussian curvature for surfaces in M 3 , and the local isometric embedding 
problem for two-dimensional Riemannian manifolds. We prove a general local existence 
result for a large class of Monge- Ampere equations in the plane, and obtain as corollaries 
the existence of regular solutions to both problems, in the case that the Gaussian curvature 
possesses a nondegenerate critical point. 
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1. Introduction 

Let K(u, v) be a function denned in a neighborhood of a point in M 2 , say 
(u, v) — 0. A well-known problem is to ask, when does there exist a piece of a surface 
z = z(u,v) in R 3 having Gaussian curvature Kl 

The classical results on this problem may be found in [10], [19], and [20]. They 
show that a solution always exists when K is analytic or K does not vanish at the 
origin. In the case that K > and is sufficiently smooth, or K(Q) = and V-fT(O) 
0, C.-S. Lin provides an affirmative answer in [15] and [16] (see [4] for a simplified proof 
of [16]). When K < and VK possesses a certain nondegeneracy, Han, Hong, and 
Lin [8] show that a solution always exists. Furthermore, if K degenerates to arbitrary 
finite order on a single smooth curve, then Q. Han and the author independently 
provide an affirmative answer in [5] and [11] (see also [6] for improved regularity). 
For an excellent survey of these results and related topics, see [7]. In this paper we 
prove the following, 

Theorem 1.1. Suppose that the origin is a nondegenerate critical point for K and 
K G C , I > 100 . Then there exists a piece of a C l ~ 9S surface in M. 3 with Gaussian 
curvature K . 

If a surface in M 3 is given by z = z(u, v) , then its Gaussian curvature is given by 

z uu z vv -zl=K(l + \Wz\ 2 ) 2 . (1.1) 
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Therefore our problem is equivalent to the local solvability of the above equation. 

Another well-known and related problem, is that of the local isometric embedding 
of surfaces into M 3 . That is, if (M 2 ,ds 2 ) is a two-dimensional Riemannian manifold, 
when can one realize this, locally, as a small piece of a surface in IR 3 ? Suppose 
that ds 2 = Edu 2 + 2Fdudv + Gdv 2 is given in the neighborhood of a point, say 
(u,v) = 0. Then we must find three function x(u,v), y(u,v), z(u,v), such that 
ds 2 = dx 2 + dy 2 + dz 2 . The following strategy was first used by J. Weingarten [25]. 
We search for a function z{u,v), with \Vz\ sufficiently small, such that ds 2 — dz 2 
is flat in a neighborhood of the origin. Suppose that such a function exists, then 
since any Riemannian manifold of zero curvature is locally isometric to Euclidean 
space (via the exponential map), there exists a smooth change of coordinates x(u, v) , 
y(u, v) such that dx 2 + dy 2 = ds 2 — dz 2 . Therefore, our problem is reduced to finding 
z{u,v) such that ds 2 — dz 2 is flat in a neighborhood of the origin. A computation 
shows that this is equivalent to the local solvability of the following equation, 

(z u - r[ lZi ) (z 22 - r 22 z t ) - ( Zl2 - r 12 zi) 2 = k(eg -f 2 - e z \ - Gz\ + 2F Zl z 2 ), (1.2) 

where Z\ = dz/du, z 2 — dz/dv, Zij are second derivatives of z, and Y l - k are Christof- 
fel symbols. For this problem we obtain a similar result to that of theorem 1.1. 

Theorem 1.2. Suppose that the origin is a nondegenerate critical point for K and 
ds 2 G C l , I > 102. Then there exists a C l ~ 100 local isometric embedding into M 3 . 

We note that A. V. Pogorelov has constructed a C 2,1 metric with no C 2 isometric 
embedding in M 3 . Other examples of metrics with low regularity not admitting a local 
isometric embedding have also been proposed by Nadirashvili and Yuan [17]. Fur- 
thermore, an alternate method for obtaining smooth examples of local nonsolvability, 
for equations with similar structure, may be found in [12]. 

Equations (1.1) and (1.2) are both two-dimensional Monge- Ampere equations. 
With the goal of treating both problems simultaneously, we will study the local solv- 
ability of the following general Monge- Ampere equation 

det(zij + a>ij(u, v, z, Vz)) = Kf(u, v, z, Vz), (1.3) 

where a,ij(u,v,p,q) and f(u,v,p,q) are smooth functions of p and q, f > 0, and 
a,ij(0,0,p,q) = d a aij(0, 0, 0, 0) = 0, for any multi-index a in the variables (u,v) 
satisfying |a| < 2. Clearly (1.1) is of the form (1.3), and (1.2) is of the form (1.3) if 
r* fe (0) = 0, which we assume without loss of generality. We will prove 

Theorem 1.3. Suppose that K(0) = \VK{0)\ = 0, det Hessi^O) ^ or HessK(O) 
has at least one negative eigenvalue, and K , a^, f G C l , l> 100. Then there exists 
a C l ~ 98 local solution of (1.3). 

Remark. 1) The methods carried out below may be slightly modified to yield the 
same result for the case when Hessi^(O) has at least one positive eigenvalue; and 
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therefore ultimately include the case of genuine second order vanishing, that is, when 
K(0) = \VK(0)\ = and |V 2 i^(0)| ^ 0. It is conjectured that local solutions exist 
whenever K vanishes to finite order and the aij vanish to an order greater than that 
ofK. 

2) Recently Q. Han [9] together with the author, have shown that local solutions 
exist for the isometric embedding problem whenever K vanishes to finite order and 
the zero set K^ifS) consists of Lipschitz curves intersecting transversely at the origin. 
Unfortunately the methods of [9] breakdown when the transversality assumption is re- 
moved. Therefore theorem 1.3 {which allows tangential intersections) and the methods 
used to prove it, may be considered as a first step towards the general conjecture. 

Equation (1.3) is elliptic if K > 0, hyperbolic if K < 0, and of mixed type if 
K changes sign in a neighborhood of the origin. Furthermore, the order to which K 
vanishes determines how (1.3) changes type in the following way. If K(0) = and 
V-fT(O) 7^ [16], then (1.3) is a nonlinear perturbation of the Tricomi equation: 

VZuu i Z vv 0. 

In our case, assuming that the origin is a nondegenerate critical point for K , (1.3) is 
a nonlinear perturbation of Gallerstedt's equation [3]: 

if Z uu i Z vv 0. 

Therefore, if sufficiently small linear perturbation terms are added to the above two 
equations, then the first (second) partial w-derivative of the z uu coefficient will not 
vanish for the Tricomi (Gallerstedt) equation. It is this fact, which allows one to 
obtain appropriate estimates for the linearized equation of (1.3) in both cases. This 
observation, lemma 2.3 below, is the key to our approach. 

From now on we only consider the case when Hessi^(O) has at least one negative 
eigenvalue, since the case of two positive eigenvalues may be treated by the results 
in [15] when K is nonnegative. Therefore, we can assume without loss of generality 
that 

Kf(u,v,z,Vz) = -v 2 + 0(\u\ 2 + \v\ 3 + \z\ 2 + \Vz\ 2 ). 

Let e be a small parameter and set u = e A x, v = e 2 y, z = u 2 /2 — v 4 /12 + e 9 w. Then 
substituting into (1.3) and cancelling e 5 on both sides, equation (1.3) becomes 

- y 2 w xx + w yy + eF(e, x, y, w, Vw, V 2 w) = 0, (1.4) 

where F(e,x,y,p,q,r) is smooth with respect to e, p, q, and r. Choose xq, yo > 
and define the rectangle X = {(x,y) \ \x\ < xq, \y\ < yo} . Let ip G C°°(X) be a 
cut-off function such that 



and cut-off the nonlinear term by F(e, x, y, w, Vw, V 2 u>) = ipF . Then solving 



= — y 2 w xx + w yy + eF(e, x, y, w, Vw, V 2 u>) = in X, 



(1.5) 



is equivalent to solving (1.3) locally at the origin. 

In the following sections, we shall study the linearization of (1.5) about some 
function w. The linearized equation is a small perturbation of Gallerstedt's equation, 
which as mentioned above admits certain estimates. These estimates are sufficient for 
the existence of weak solutions, however the perturbation terms cause some difficulty 
in proving higher regularity. To avoid this problem, we will regularize the equation 
by appending a suitably small fourth order operator. In section §2 we shall prove the 
existence of weak solutions for a boundary value problem associated to this modified 
linearized equation. Regularity will be obtained in section §3 . In section §4 we make 
the appropriate estimates in preparation for the Nash-Moser iteration procedure. 
Finally, in §5 we apply a modified version of the Nash-Moser procedure and obtain 
a solution of (1.5). 

2. Linear Existence Theory 

In this section we will prove the existence of weak solutions for a small perturbation 
of the linearized equation for (1.5). Fix a constant A > 0, and for all i,j = 1,2 let 
bij, h, be C r (R 2 ) be such that: 

i) the supports of b^ , bi, and b are contained in X, and 

w) J2 \bij\c 10 + \bi\c 10 + \b\c 10 < A. 
We will study the following generalization of the linearization for (1.5), 



where x l = x, x 2 = y and a u = -y 2 + eb u , a 12 = eb 12 , a 22 = 1 + eb 22 , ai = eb x , 
a 2 = eb 2 , a = eb. 

To simplify (2.1), we shall make a change of variables that will eliminate the mixed 
second derivative term. In constructing this change of variables we will make use of 
the following lemma from ordinary differential equations. 

Lemma 2.1 [1]. Let G(x,t) be a smooth real valued function in the closed rectangle 
\x — s\ <T\, \t\ < T 2 . Let M = sup \G(x,t)\ in this domain. Then the initial-value 
problem dx/dt = G(x,t), x(0) = s, has a unique smooth solution defined on the 
interval \t\ < min(T 2 ,Ti/M) . 

We now construct the desired change of variables. 




(2.1) 
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Lemma 2.2. For e sufficiently small, there exists a C r diffeomorphism 

€ = t(x,y),ri = y 
of X onto itself, such that in the new variables (£, rj) 

L = y^"ij<lr..r, +^2~aid Xi + a, 

where_x 1 = x 2 = r]_, a u = —r} 2 + £ bu, «i2 = 0, a 2 2 = 1 + ^22, o~i = £b~i, a 2 = eb 2 , 
a — eb, and bij , h, b satisfy: 

i) b ij ,b i ,beC r - 2 (X) ) 

ii) bij , b~i , and b vanish in a neighborhood of the lines £ = ±x , and 

in) E \bij\cs(x) + \bi\cHx) + \b\c»(x) < A '> 
for some fixed A' . 

Proof. Using the chain rule we find that 012 = 0126c + 022^ ■ Therefore, we seek a 
smooth function £(x,y) such that 

ai2^ + a 2 2^ = in X, £(x,0) = x, £(±x , y) = ±x . (2.2) 

The boundary condition ^(±xo,y) = ±xq states that the vertical sides of dX will be 
mapped identically onto themselves under the transformation (£, 77) . Moreover, the 
horizontal portion of dX will be mapped identically onto itself since rj — y. Thus, 
(C,rj) will act as the identity map on dX . 

Since 012 = £612 and a 2 2 = 1 + £b 22 , by property (ii) if e is sufficiently small 
the line y = will be non-characteristic for (2.2). Then by the theory of first order 
partial differential equations, (2.2) is reduced to the following system of first order 
ODE: 



X = 


Ol2 

i 


x(0) 


= s, 


— Xq < S < X , 
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y = 


1, 


1/(0) 


= 0, 




i = 


0, 


m 


= s, 


£(±xo,y) = ±xo 



where x = x(t), y = y(t), £(t) = £(x(t),y(t)) and x, y, £ are derivatives with 
respect to t. 

We first show that the characteristic curves, given parametrically by (x,y) = 
(x(t),t), exist globally for — yo < t < y$. We apply lemma 2.1 with T\ = 2xq and 
T 2 = yo to the initial- value problem x — — , x(0) = s. By property (ii) for the by 



so for e small, M < Thus min(T 2 , Ti/M) = yo, and lemma 2.1 gives the desired 
global existence. 

We observe that £ = s is constant along each characteristic. In particular, since 
mi(±a; 0)2/ ) = the characteristics passing through (±x ,0) are the vertical lines 
(±xo,t), so that £(±xo, y) = is satisfied. 

We now show that the map p : X — > X given by 

(s,t) I— >■ (x(s,t),?/(M)) = 

is a diffeomorphism, from which we will conclude that £ = s(x, y) is a smooth function 
of (x,y). To show that p is 1-1, suppose that p(si,£i) = p(s 2 ,£ 2 ). Then t x = t 2 and 
— ^(52,^2), which implies that Si = s 2 by uniqueness for the initial-value 
problem for ordinary differential equations. To show that p is onto, take an arbitrary 
point (x\,yi) G X, then we will show that there exists s G [— Xq,Xo] such that 
p(s,yi) = (x(s,yi),yi) = (x 1 ,y 1 ). Since the map x(s r ) : [-x ,x ] ->■ [-x ,a:o] is 
continuous and x(±xo, •) = ±^o ? the intermediate value theorem guarantees that 
there is s G [— x ,x ] with = xi, showing that p is onto. Therefore, p has a 

well-defined inverse. 

To show that is smooth it is sufficient, by the inverse function theorem, to 
show that the Jacobian of p does not vanish at each point of X . Since 



Dp 



X s Xf 

o 1 



this is equivalent to showing that x s does not vanish in X . Differentiate the equation 
for x with respect to s to obtain, -j t {x s ) = (^) x x S) x s (0) = 1. Then by the mean 
value theorem 



\x s (s,t) - 1| = \x s (s,t) - x s (s,0)\ < y sup \(— ) x \ sup \x, 

X 0,22 X 

for all (s,t) G X. Thus, by property (ii) for the 6^, 

1 - £Ciy sup \x a \ < x s (s,t) < eCiyosup \x s \ + 1 
x x 

for all (s,t) G X. Hence for e sufficiently small, x s (s,t) > in X. We have now 
shown that p is a diffeomorphism. Moreover, by lemma 2.1 and the inverse function 
theorem we have p, G C r . 

Lastly, we calculate a n , a 22 , ai, a 2 , and show that they possess the desired 
properties. It will first be necessary to estimate the derivatives of £ . By differentiating 
(2.2) with respect to i, we obtain 

(—)&)* + (e.)» = -(— u*, o) = i. 

o 22 a 22 
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As above, let (x(t),y(t)) be the parameterization of an arbitrary characteristic, then 
£x00 = €x(x(t),y(t)) satisfies f x = -(^) x f x , f x (0) = 1. By the mean value theorem 



k 0122 

- 1| = \Ut) - U0)\ < yosup |(— ) x | sup \u 

X 0>22 X 

By property (ii) for the by , 

1 - eCxyoSup |f x | < £ x (t) < eCiy snp |f x | + 1. 

X X 

Since this holds for any characteristic, we obtain 

sup |f x | < l — — := C 2 . 

x 1 - edyo 

It follows from (2.2) that 

sup If j, I < C 3 , 

X 

where C 2 , C3 are independent of e and hij. In order to estimate f xx , differentiate 
(2.2) two times with respect to x: 

(^)(e«)x + (ixx)y = ~2{^U XX - (^)Ux, U(X,0) = 0. 
022 °22 a 22 

Then the same procedure as above yields 

sup |f xx | < edyo sup |f xx | + eC 5 y , 
x x 

implying that 

sup f xx < — := eC 6 . 

x 1 - eC A y 

Furthermore, using the above estimates we can differentiate (2.2) to obtain 

sup < eC 7 , sup \£ m \ < eC 8 , 

X X 

for some constants Cj , C$ independent of e and b^ . This procedure may be continued 
to yield, 

\d a t\ < eC 9 , 

for any multi-index a satisfying 2 < \a\ < 10. 

We now show that an , a 2 2 , a±, a 2 satisfy properties (i) , (ii) , (Hi) and have the 
desired form. Calculation shows that, 

an = a ll£ x + 2ai2fa;£y + a 22^, °1 = a ll£,xx + 2ai2^a;j/ + a 22^yy + Oif x + ^2^- 
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Furthermore, according to the above estimates and the fact that the bij vanish in a 
neighborhood of dX , we may write 

ix = 1 + ex, 

where x £ C r_1 (X) vanishes in a neighborhood of the lines x = ±x . It follows that, 

a n = -if + eb u , ai = ebi, 

where b~u and b~i satisfy properties (i) , (ii) , (iii) . Moreover, since 022 = 022 and 
a 2 = a 2 , properties (i) , (ii) , (Hi) hold for these coefficients as well. □ 

For the remainder of this section and section 3, (£, 77) will be the coordinates of 
the plane. For simplicity of notation we put x — £, y — rj, and = a^, a« = a,, 
a = a , % = % , k = bi, b = b. 

In order to obtain a well-posed boundary value problem, we will study a regular- 
ization of L in the infinite strip fl = {(x,y) \ \x\ < x } . More precisely, define the 
operator 

Lg = —Od XX yy + L, 

where > is a small constant that will tend to zero in the Nash-Moser iteration 
procedure. Furthermore, we will need to modify some of the coefficients of L away 
from X as follows. First cut b^ , bi , and b off near the lines y = ±yo , so that by 
property (ii) of lemma 2.2 these functions vanish in a neighborhood of dX , and the 
coefficients a^-, a^, and a are now defined on all of Q. Choose values y±, yi-, and 
y 3 such that yo < y\ < yi < y 3 , and let 5 > be a small constant that depends on 
2/2 — l/i and y 3 — y 2 . Then redefine the coefficient a in the domain Q — X so that: 

i) ae C r - 2 (n), 

ii) a = 1 if 1 2/ j > yi, 

iii) a > for \y\ > y , 

iv) d y a > if y > y , and d y a < if y < —yo- 
Redefine an in f2 — X and near <9f2 so that: 

i) an G C r - 2 (H), 



ii) an 



-y 2 if l/o < \y\ < Vi, 

_(lfi±ia)2 if|y|>y 2 , 



iii) (9 y aii < if y > yo, and d y au > if y < — yo, 
if) sup^ dyydu < 5 , 
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v) an Ian < -d, d%a>u\m = 0, a <r - 2, and sup n |dfan| < eA' , 1 < /3 < 8. 
Lastly, redefine a 2 in O - X so that: 
i) a 2 G C r - 2 (U), 

if yo < 1 2/ 1 < 2/2, 

U) a 2 = { -Sy + if 2/ > 2/3, 



if < -2/3, 



m) a 2 < if 2/ > 2/2 ! and a 2 > if 2/ < —2/2 , 
iv) sup| w |> w |9 v a 2 | < S. 
Denote the operator L with coefficients modified as above by L' , and define 

Lq = —9d xxyy + L'. 

Note that since we are studying a local problem, as stated in the introduction, we 
may modify the coefficients of the linearization away from a fixed neighborhood of 
the origin. This will become clear in the final section, where a modified version of 
the Nash-Moser iteration scheme is used. 

Consider the following boundary value problems 

L e u = f in Q, u\ dQ = 0, (2.3) 

L d u = f in Q, u x \ d n = 0, (2.4) 
and the corresponding adjoint problems 

L* 9 v = g in Q, v\ dn = 0, (2.5) 

L* e v = g in Q, v x \ dn = 0, (2.6) 

where L* e is the formal adjoint of L$. The main result of this section is to obtain 
weak solutions for all four problems. 

We will make extensive use of the following function spaces. For m,n G Z> let 

C (m ' n) (Q) = {u : n ->■ R | fi£fi£u G C°(n), a < m, (3 < n}, 
& m ' n) {tt) = {u G C (m ' n) (n) | u\ m = 0, u has bounded support}, 
^'"'(H) = {uG C (m ' n) (H) | u^lan = 0, u has bounded support}. 
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Define the norm 

II « H (m> „)= ( E II ^« iii»(n)) 1/2 , 

a<m, /3<n 

and let H^ m ' n \ft) and H { x m ' n) (VL) be the respective closures of C (m ' n) (Tl) and C x m,n) (9) 
in the norm || • ||( m ,n)- Furthermore, let H m (Vt) denote the Sobolev space of square 
integrable derivatives up to and including order m, with norm || • || m . Denote the 
L 2 (f2) inner product and norm by (•, •) and || • || respectively, and define the negative 
norm 

II II 1(^)1 

II U ||(_ m _ n )= sup 1| 1| . 

vEH( m > n ) (fi) II V \\(m,n) 

Let H^ m ^ n \Vt) be the closure of L 2 {VL) in the norm || • || ( - m - n ), then H^-^iVL) 
is the dual space of H^ m ' n ^ (Q) . The dual space of H x m ' n ^ (Q) is defined similarly. 

Let / G L 2 (VL). A function u G L 2 (Q) is said to be a weak solution of (2.3) 
(respectively (2.4)) if 

(u,L* g v) = (f,v), for all v G C°°(Il) (for all v G C™(Ti)). 

We shall employ the energy integral method, developed by K. O. Friedrichs and 
others, to prove the existence of weak solutions for (2.3) and (2.4). The first step is 
to establish an a priori estimate. 

Lemma 2.3 (Basic Estimate). If e , 9 , and 5 are sufficiently small, then there exist 
constants C\,C2 > independent of e , 6, 5, and functions A, B,C, D, E G C°°(f2) 
where E > and E = 0(\y\) as \y\ — > oo , such that: 

(Au + Bu x + Cu y + Duyy, L gu) > 

Cl[\\ U || 2 + || EUy || 2 +9(\\ U X || 2 + || U X y || 2 + || Uyy ^ + 6 \\ U X yy || 2 )], 

for all u G C°°(f2) with bounded support such that u x (—x ,y) = 0, and either 
u(xo,y) = or u x (xo,y) = 0. Furthermore, 

|| U || + || U„ || +V0(\\ ^x || "I - || ^" X y || "I - || ^"yy || + V / # || ||) < C 2 || L e M ||, 

for all u G C^itt) and for all u G C™(tt). 

Proof. We first define the functions A,B,C and D. Let /i be a positive constant 
such that \n + a±i > 1 throughout fi, and let 7 G C°°([— rr , ^0]) be such that 

7 (x) = 



1 if -Xq < X < ^, 

ifx = x"„, " " 



with 7(2;) > except at x = x , and 7' < 0. Define 
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c = 



fxd y a n if \y\ < y , 
-2/iy if \y\ > y , 

D = 9, 



and note that A, B,C,D e C°°(fi) . 

We now prove the first estimate. Let u G C°°(fi) satisfy the given hypotheses. Let 
(711,712) denote the unit outward normal to dVL. Then integrate by parts to obtain: 

(Au + Bu x + Cu y + Duyy, L e u) = 




hU 2 x yy + h^yy H yyU ^ ^" ^4 U 3 y ^" X yU ^ 

+ 2I 6 U X yUy + I 7 UI + 2IgU x Uy + IfjUy + 1^ 

+ / -A^xj/ + -hu xy u x + J 3 m^ + J 4 m^ + J 5 m 2 ; 
an 



where 



Ji = -OBri!, J 2 = OByUi, J 3 = ^-BonTii, 

2J 4 = — 9A x rii — 9C xy rii + (Dau) x ni — Daini, 

2J 5 = —(Aau) x ni + Aa\n\ + Ban\ + 9A xyy n\, 

and the remaining Ji, . . . , J 10 will be given below as each term is estimated. First 
note that J 2 \on = J^dn = 0. Furthermore J\ = ■ ■ ■ = J 5 = on the portion of 
the boundary, x = Xq, since 7(^0) = 0. Whereas on the other half of the boundary, 
x = —xq, we have u x (— xo,y) = and J5 = \Bari\ > 0. It follows that the entire 
boundary integral is nonnegative. 

We now proceed to estimate the integral over Q , beginning with 7i , I 5 , and I w , 
which are given by 

h = 6D, J 5 = —-9B y , 

2/io = (Aaii) xa! + {Aa 22 )yy - (Aa^x - {Aa 2 ) y 

+ 2Aa - (Ca) y - {Ba) x - 6A xxyy + (Da) yy . 

Since B is a function of x alone, 7 5 = , and by definition of D , 1\ = 6 2 . It will now 
be shown that I w > M x in Q, for some constant Mi > independent of e and 9. 
In order to accomplish this we shall treat the regions \y\ < y , y < \y\ < yi, 
Hi < \y\ < U2, and \y\ > y 2 separately. Moreover, throughout this proof Mj, 
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i = 1,2,..., will always denote positive constants independent of e and 9. A com- 
putation yields, 

1 1 

ho = -a 22 d yy a u - a u a - -Cd y a - -(Aa 2 ) y + 0(e + 6). 

In the region \y\ < y we have a,d y a,a2,d y a 2 = 0(e), a 22 = 1 + 0(e), and 
d yy CLii = — 2 + 0(e), so that here ho > M 2 . If yo < |y| < yi, the conditions 
placed on a guarantee that 

—ana — -Cd y a > 0; 

furthermore a 22 , an , and a 2 have the same properties in this region as in the previous. 
Hence, I w > M 3 when y < \y\ <y x . If yi < \y\ < y 2 then 

-a 22 d yy a n = 0(5), -a n a>yj, a 2 = d y a = 0, 

showing that ho > M 4 in this region. Lastly, when \y\ > y 2 we have ho > M 5 since 

d yy a n = d y a = 0, - a u a = ( Vl+ 2 V2 ) 2 , - ^(Aa 2 ) y = 0(5). 

The desired conclusion now follows by combining the above estimates. 
Next we show that 

/ / ^ u m + 2 hu yy u xy + hu 2 xy > M 6 6(\\ u yy f + || u xy f), 
J Jn 

where ^ 

h = ~-jQD xx + Da 2 2, ^3 = --j^C x , 



1 = -\ec y - l -eB x -6A + D ail . 



This will follow if I 2 > M 7 9 , J 4 > M 8 9 , and I 2 h — If > . A calculation shows that 

I 2 = 6a 22 = 9(1 + 0(e)), h = 0(e9), 

h = 29(a u - hj y ) + 0(e9) = 29(fi + a n + 0(e)). 

Therefore, since was chosen so that fi + an > 1 in Q, the desired conclusion follows 
if e is sufficiently small. 
We now show that 



I In 



I 7 u 2 x + 2I s u x u y + I 9 u 2 y > Mg(9 || u x \\ 2 + || Eu y || 2 ), 

n 

where 

2I 7 = -2Aa n - (Ba n ) x + 25a! + (Ca n ) y + 9B xyy + 9A yy - (Da n ) yy , 
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27 8 = -(Ba 22 ) y + Ba 2 - (Ca u ) x + Ca x + 9A xy + (Da u ) xy - (Da^y, 

27 9 = -2Aa 22 - (Ca 22 ) y + 2Ca 2 + 6C xxy + 6A XX 

- (Da u ) xx - {Da 2 ) y + (Dai) x - 2Da. 

Again, this will follow if I 7 > M 1Q 0, I 9 > M n E 2 , and I 7 I 9 - If > 0. A calculation 
shows that 

I i 

h = a\i + 2 Cd v an + °(- d yy a n + ipx^u + 0(e)), 

11 11 

1 8 = --C x a u ~ ~Cd x a u + -Cai + -Ba 2 + 0(9), 

1 9 = (a u - C y )a 22 + Ca 2 + 0(e + 9) 

= (2/i + a n + 0(e))(l + 0(e)) + Ca 2 + 0(e + 9). 

Then 7 9 > M U E 2 immediately follows since Ca 2 = 0(e) if \y\ < y , Ca 2 > if 
\y\ > Vo, Ca 2 = 0(\y\ 2 ) as \y\ — > oo, and 2/i + an > I. To show that J 7 > M w 9, we 
consider the regions \y\ < yo and \y\ > yo separately. If \y\ < yo then 

CdyCiu = fi(d y a u ) 2 > 0, - 9 sy oii = 2 + 0(e), 7 x a n > -0(e), 

so that here I 7 > 29 + 0(e9). Furthermore, when |y| > y we have J 7 > y\ + 0(0) 
since 

a?i > Vo, Cd y au > 0. 
Finally, I 7 Ig — i| > follows from the next calculation. If \y\ < y then 

hh ~ I\ > (a 2 u + ^(d y a n ) 2 + 29 + 0(e9))(l + 0(e + 9)) 

- \0(e 2 )a 2 u - ±0(e 2 )(d y a u ) 2 - 0(e9 + 9 2 ), 

whereas if \y\ > y then 

hh - if > (vt + 0(9))(1 + 0(5y 2 )) - 0(9 2 y 2 ). 

Lastly, we deal with the term 2I 6 u xy u y . Consider the quadratic form: 

M 6 9u 2 xy + 2I 6 u xy u y + M 9 E 2 u 2 y , 

where ^ 

h = ~^Ba 22 . 

Since 

(M 2 9)(M 3 E 2 ) - I 2 > M n 9 - M 12 9 2 (l + 0(e)) 
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for some M n , M 12 , we obtain 

M 6 6u 2 xy + 2I e u xy u y + M 9 E 2 u 2 y > M l3 (6u 2 xy + E 2 u 2 y ). 

This completes the proof of the first estimate. 

To obtain the second estimate we need only observe that the above arguments hold 
if B = and u G C°°(Q) or u G C^fi). Then an application of Cauchy's inequality 
(ab < \a 2 + j^b 2 , A > 0) yields the desired result. The reason for including B in the 
first estimate will soon become clear. □ 

Having established the basic estimate, our goal shall now be to establish dual 
inequalities of the form: 

II v || < Ci || L* e v || (-1,-2) for all v G C°°(n), 

|| v || < C 2 || LJu II (-i-2) for all v G 

The existence of weak solutions to problems (2.3) and (2.4) will then easily follow 
from these two dual estimates, respectively. In order to establish the dual estimates, 
we will need the following lemma. Let P denote the differential operator 

P = Dd 2 y + Bd x + Cd y + A, 

where A, B, C, and D are defined in lemma 2.3. Note that P is parabolic in f2, away 
from the portion of the boundary, x = x . This is the reason for including B in the 
first estimate of lemma 2.3. 

Lemma 2.4. For every v G C°°(Q) there exists a unique solution u G C°°(fl) D 
H 4 (tt) C C°°(fi) n C 2 (tt) of 

Pu = v in Q, u(-x ,y) = u x (-x ,y) = 0, u(z o ,y) = 0. 

Furthermore, for every v G C£°(f2) £/iere exzsfo a unique solution u G C°°(f2) D 
H\n) C C°°(fi) nC 2 (H) of 

Pu = v in Q, u x (-x ,y) = 0, w^o, y) = 0. 

Proof Let r > be a small parameter, and define the subdomains 

f2 r = {(x, y) | -x < x < x - t}. 

Then P is parabolic in fl T for each r. 

We now consider the case when v G C°°(Q). The parabolicity of P guarantees 
the existence (see [13]) of a unique solution to the Cauchy problem 

Pu = v in Q, u(—x ,y)=0, 
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such that u G H°°(Q T ) for every r. Furthermore, u x (— x ,y) = since 

Bu x \(- XOi y) = Pu\^ x , hy) = v(-x ,y) = 0. 

We shall now show that u G H A (Q) . This will be accomplished by estimating the 
H 4 (Q T ) norm of u in terms of the H 4 (Q) norm of v, independent of r. To facilitate 
the estimates, we first construct an appropriate approximating sequence, {u k } ( k x L 1 , 
for u. Define functions u k G C°°(R) by 





if 


\y\ 




if 


\y\ 



Mv) = < n J" : n ; (2.7) 



such that < v k < 1, sup \v' k \ < |, and \vk\c*(n.) < M for some constant M 
independent of k. Let u k = v k u, then: 

i) u k G C°°(U T ) for all r, 

n) u k has bounded support and u k (—x ,y) = u k (—x ,y) = 0, 

Hi) \\ u — u k ||4,n T — * as k — > oo, 

if) || C-Uj, — Cw^ ||n T — >■ as fc — >■ oo, 

where C was defined in lemma 2.3. All of the above properties are evident except for 
(iv) , and (iv) follows from the following calculation. Let 

ni k ^ = {(x,y)en T \k 1 <\y\<k 2 }, 

then 

ll CUy ~ CU k \\l T < || C(Uy ~ V k Uy) f + || CVfc 



< 



I I ° 2 < + I I (« 2 (r) 2 ^ 



where \x was defined in the proof of lemma 2.3. By solving for Cu y in the equation 
Pu = v , we have 

Cu y = v — Duyy — Bu x — Au G L 2 (f2 r )- 

Therefore 

C 2 ul ->■ as fc ->■ oo. 

a (fc,oo) f 
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Furthermore, 

36/i 2 -u 2 — > as k — > oo 




since u G L 2 (Q T ). This proves (iv) . 

We now proceed to estimate the H 4 (il T ) norm of u. Let ( = ((y) G C°°(R) be 
such that C < 0, C(y) = -\y\~ 1/2 if \y\ > Vu C'(y) > if y > 0, and c'(y) < if 
y < 0. Then set k = 2 sup | C a n I > an d integrate by parts to obtain 

(Ku k yy + (u k )Pu k =jj [ K D}(u k yy ) 2 + [-D( + ^B x -^C y -A)](u k ) 2 

+ [\^ yy + \{DQ yy - 1 -{BQ X - l -{CQ y + (A](u k ) 2 

The boundary integral is nonnegative since u k (—x ,y) = u k (—x ,y) = 0, and 
-KSni|( a . _ T)W ),SCni|( X0 _ T)1/ ) > 0. Also, kD > 0, 




--DC + <\b x - ±C„ - A) > k(2h + a n + 0(e + 6)) > k, 



and 



\**vu + \(DQ m - 1 -{BQ X - l -{CQ y + (A 

11 11 

= -^ Kd yv a u ~ if^Cy - C^n + ^{PQ m - -{BQ X + 0(e) 

{K-( ail + 0{e + 6) if|y|<yi, 
|2/I~ 1/2 [|M + an + 0(6)] + 0(k5) if y x < \y\ < y 2 , 
\ y \-V*[lp + ail + 0(e)] if \y\>y 2 . 

Therefore if e, 9, and 5 are sufficiently small, we may apply the Schwartz inequality 
followed by Cauchy's inequality to obtain 



-(u* \\n T + || u; \\n T + || u K yy ||n T < Mi || Pu" \\n T , 

for some constant Mi independent of r. The properties of u k guarantee that by 
letting k — > oo, we obtain 

(u \\q t + || u y ||n T + || u yy \\ Qt < M x \\ Pu ||n T = Mi || v \\q t < M x \ v\. 

We now estimate d x d^u for a = 1, . . . , 4, and j3 = 0, 1, 2. Differentiate Pu = v 
with respect to x: 

D(u x ) yy + B(u x ) x + C(u x ) y + (A + £ x )u x = f x - Cx^ - A x u. (2.8) 
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Since u x (— x ,y) = and A x , C x vanish outside a compact set, we can apply the 
same procedure as above to obtain 



~C u x \\n T + || u xy ||n T + || u xyy ||n T < Mi || v x - C x Uj, - A x u \\ Qt 

< M 2 (\\ v x \\n T + || u y \\n T + || u \\n T ) 

< M 3 (|MI + IK ID- 
Differentiating (2.8) with respect to x produces 

D{u xx )yy + B(u xx ) x + C{u xx )y + (A + 2B x )u xx 

= v xx ~ d x (C x u y + A^m) — C x u xy — [A x + B xx )u x := fi. 

Again we apply the same method. However, since u xx {— xo,y) = B~ 1 v x \(- XOtV ) from 
(2.8), we now have 

|| a/~C M ^ lbr + ll U xxy ||fi T + || Uxxj/j/ ||fi T < M || ^ ||f} T +M 4 

< M 5 (|| u ll + H ^ ll + H ||) + M 4 , 

where M 4 = n\B\~ 1 (f x= _ xo v xy + f 2 ) 1 / 2 which is independent of r. We can estimate 
II V~(dx u lln T ) « = 3,4, and || d x d!^u \\q t , a = 3, 4, /3 = 1, 2, in a similar manner. 
To estimate , differentiate Pw = v with respect to y : 

D{Uy)yy + B{Uy) x + C{Uy)y + (A + Cy)Uy =Vy~ AyU . (2.9) 

Since u y (—x ,y) = 0, C y < 0, and A, vanishes outside a compact set, the same 
method as above yields 

|| V~CUy \\n T + || Uyy \ | ^ + || Uj,^ | | ^ ^ Ml || Vy ~ AyU \ | ^ 

< M 6 (||^|| + || 

Furthermore, || tta^y ||n T an d || ""^^ ||q t can be estimated by differentiating (2.9) 
with respect to x and y, respectively. 

The combination of all the above estimates produces, 



£ II V^Cdy |k + £ II |k< M 7 || „ || 4 +M 8 , 



O=0 a + /3<4 



where M7 and Mg are independent of r. Then letting r — > we find that 
d%d$u G L 2 (fi), a + /3 < 4, /3 ^ 0, and that yf^d^u e L 2 (ty, a = 0,...,4. 
It follows that u G H 4 (K) for every compact if C f2, so that u G C 2 (Q). 

We now show that d x u G L 2 (fi) , ct = 0, . . . , 4. Let Qi, q 2 G C 00 (R) be given by 

, x + if -x < x < =f- , f-y if < 2/0, 

10 if < x < xq. 10 if \y\ > T, 
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such that £2(2/) < if y > and £2(2/) > if y < 0, where T > is large enough so 
that — 1 < g 2 < e. Then define B — B + Qi and C = C + £2 — Sfj,d y bu = — + £2 , 
and set 

P = Bd x + Cd y + A. 
If w G C£°(Q), then integrating by parts yields 

{w,Tw) = I iH^ ~ & + A]w2 + L [ ~^ ni]w2 ' 

The boundary integral is nonnegative since B(—xo,y) = 9 and B(x ,y) = 0. Fur- 
thermore, 

~B X -^C y + A = -q 2 - a u + 0(e + 9) > M 9 

for some constant M 9 > 0. Thus 

II w ||< M10 || P*w || . (2.10) 

Since i> — -Dm^ + + (02 — £^d y bn)u y G L 2 (Q), (2.10) implies (see the proof of 
theorem 2.1 below) the existence of a weak solution u G L 2 (Q) of 

Pu = v - Du, yy + Qiu x + (g 2 - efxd y b u )uy, u(-x , y) = 0. 

We shall now show that u = u. Since P is a first order differential operator, 
we may apply G. Peyser's extension [21] of Friedrichs' result [2] on the identity of 
weak and strong solutions to obtain a sequence {u k }'j? =1 , such that u h G C°°(f2) has 
bounded support, satisfies u k (—x ,y) = 0, and 

|| u — u k || + || Pu k — (v — Dtiyy + Q\U X + (g 2 — £fidybii)u y ) \ \ — > as k — > oo. 

Set v k = u — u k . Using the fact that |y|~ 1 / 4 f fe — > |y| _1 / 4 (-u — u) G L 2 (f2) and recalling 
the definition of P, we have 

|(_| y |-i/V,iV)| < || |y| _1/ V IHI Pv k || 

< Mn || v - Duyy + Qiu x + (g 2 - efxdybu)u y - Pu k ||-> 0. 

Then the following calculation shows that || u — u k ||l 2 (x)— ^ for every compact 

Ken-. 



(-\y\-^v k ,Pv k ) = lim / / [Jlyl-^ + idyl-V^-lyr 1 /^]^) 2 

+ I [-hy^Cn.-^y^Bn^f 
> lim f I l\y\- 1/ \\ f i + a 11 - 1 - + 0(e + 9))](v k y 



> M 12 || \y\-'/% k \\l 



K ■ 
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Therefore u = u in L 2 (fi). 

Differentiating the equation P-u = v with respect to <9" , a = 1, . . . , 4, and apply- 
ing the above procedure, shows that d%u G L 2 (£l), a — 1, ... ,4. We now have that 
ueH 4 (tt). 

To complete the case when i> G C°°(Q), we must show that u(x ,y) = 0. Since 
B(xo,y) = 0, from the equation Pu = v we find that 

(DUyy + CUy + Au)\^y) = v(x ,y) = 0. 

Furthermore since u G H 4 (Q) , u — > as \y\ — > oo. Therefore, applying the maximum 
principle to the above equation, we have u(xo,y) = 0. 

We now consider the case when v G C™(Q). Let h(y) G if°°(R) be the unique 
solution of the ODE: 

D{-x Q , y)h" + C{-x Q , y)ti + A(-x Q ,y)h = v(-x , y). 

Then as before, the parabolicity of P guarantees the existence of a unique solution 
to the Cauchy problem 

Pu = v in Q, u(—x ,y) — h(y), 
such that u G H°°{Vt T ) for every r. Furthermore, u x (—x ,y) = since 

BU X \(- XQi y) = V(-X ,y) - {DUyy + CUy + Au) | (_ X0 t y) = 0. 

Moreover, the same methods used above can be used here to show that u G H 4 (Q). 
Lastly, to show that u x (x ,y) = 0, differentiate Pu = v with respect to x and use 
that B(xo,y) = to obtain 

(D(u x ) yy + C(u x ) y + (A + B x ) U x ) \ {xq ty ) = V x (x ,y) - (C x u y + A x u)\( XOt y) = 0. 

Since u x — > as \y\ — > oo, by the maximum principle u x (xo,y) = 0. □ 

With lemma 2.4 we are now in a position to establish the dual inequalities. 
Proposition 2.1. There exist constants M±,M 2 such that: 

|| v \\< Mi || LgV ||(-i-2) for all v E C°°(Ti), 

|| v || < M 2 || L> || (_i-2) /or a// v G 

Proo/. We first consider the case when v G C°°(n). Let u G n P 4 (fi) be the 

unique solution of 

Pu = v in n, u(-x ,y) = u x (-x ,y) = 0, u(z o ,y) = 0, 
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given by lemma 2.4. 
We now show that 

(Au + Bu x + Cu y + Duyy, L e u) > 

d[|| u || 2 + || Eu y || 2 +0(|| u x || 2 + || u xy || 2 + || m to || 2 +6 || Uaw || 2 )], 

where A, B,C, D, E , and Ci were given in lemma 2.3. Let v k be given by (2.7) and 
define the sequence {u k } < k x L 1 , where u k = v k u. Then as in the proof of lemma 2.4 we 
have: 

i) u k e C°°{Q)nH 4 {Q), 

ii) u k has bounded support and u k (—x ,y) = 0, u k (x ,y) = 0, 
Hi) || u — u k || 4 — > as k — > oo, 

if) || Eu y — Eu k ||— >■ as A; — >■ oo. 
Let {-Ufe}feli a C°° approximation of {■u fc }^ 1 such that: 

i) u k eC°°(U), 

ii) u k has bounded support and (u k ) x (— xo, y) = 0, u k (xo,y) = 0, 
Hi) || — u k 114—7- as /c — )> oo, 

if) || Eu k — E{uk)y || — )■ as A; — )■ oo. 
Then applying lemma 2.3 we have, 

(Au + -Bf X + CUy + DUyy, LqU) 

= lim (Aw fc + B(u k ) x + C(w fe ) y + D(u k ) yy , L e u k ) 

k— >oo 

> lim dHI f fc || 2 +|| E{U k )y \\ 2 +9{\\ («fc)*l| 2 +ll ("*)*» || 2 +|| ( U k)yyW 2 +0 || (U k ) X yy\\ 2 )} 

K— s-oo 

= Cl[|| U || 2 + || -Efj, || 2 +6I(|| 11, || 2 + || U X y || 2 + || Uyy ^ + ^ || U X yy \^)\. 

By the above estimate and definition of the negative norms, it follows that 

II L* e v || (_i_2) || u ||(i,2) > (L* e v,u) 

= (v,L e u) 

= (Au + Bu x + Cu y + Du y y, L gu) 
> dOlu || 2 + || || 2 

+6>(|| n x || 2 + || u XJ/ || 2 + || u w || 2 +9 || u xw || 2 )]. 
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Furthermore, using Cauchy's inequality and the equation Pu = v, we obtain 

|| L* g V ||(-l-2) > [|| U || + || EUy || +V^(|| U X || + || U X y || + || Uyy \\ + \\ U X yy ||)] 

> M ^ II v ||, 

for some constants C[,Mi > 0. Moreover, similar arguments may be used to treat 
the case when v G C£°(J1). □ 

The existence of weak solutions to problems (2.3) and (2.4) immediately follows 
from proposition 2.1 by a standard functional analytic argument. We include the 
proof here for convenience. 

Theorem 2.1. For each f G L 2 (Q) there exists a weak solution u G H^ 1,2 \Q) , 
H^ 2 \n) of (2.3), (2.4) respectively. 

Proof. We shall first treat problem (2.3). Let W = Lq(C°°(Q,)) and define the linear 
functional F : W -> R by 

F(L* e v) = (f,v). 

Using proposition 2.1, the following calculation will show that F is bounded as a 
linear functional on the subspace W of f/"( -1,-2 )(fi) , 

\F{L»\ = |(i»| 

< Mi || / mi ||(-i-2) ■ 

Use the Hahn-Banach theorem to extend F from W onto the whole space iJ( -1 '~ 2 )(fi) . 
It follows from the Riesz representation theorem that there exists u G -f/^ 1 ' 2 ^) such 
that 

F(w) = (u,w) for all w G // (_1 '- 2) (Q). 
Thus, restricting w to If we have 

(u, L* e v) = F(L* e v) = (/, v) for all veC°°(tt). 

The case of problem (2.4) may be treated in a similar manner. □ 

We now prove the existence of weak solutions for the adjoint problems (2.5) and 
(2.6). The existence of solutions for these problems will be needed in the next section, 
where they will aid in proving higher regularity for solutions of (2.3). 

The formal adjoint of Lg is given by 

L*g = —0d xxyy + Qjl\d xx + 022dyy 

+ (2d x a n - ai)d x + (2d y a 2 2 - a 2 )d y + (d xx a u + d yy a 2 2 - d x a x - d y a 2 + a). 



21 



All the coefficients of Lg, denoted a*j,a*,a*, have the same properties as the coeffi- 
cients of Lq, except a 2 = 2d y a 2 2 — a-2- This difference will not allow us to directly 
apply the above procedure to obtain weak solutions for (2.5) and (2.6). However, if 

h(x,y) = e Jo "2 2 (*.t) , 
then by setting v = hw, the equation L* e v = g becomes L* 9 w = g/h, where 

Lq @dxxyy '^@~^~&xxy 26 ^ ~j^~&xyy 



+Ki - 0^f)d xx - A6^fd xy + (a* 22 - Q h -f)d yy 

hi/ _ „ h'T.'r.'}! \ ^ / a. . ^ a. ^ T ^ n ^ 



+(a* 2 + 2a* 2 ^ - 29^)d y + (a{ + 2a* u '-f - 29^)d x 

+(aj 1 ^ + + + + a* - 6^). 

The special choice of h guarantees that the coefficient of d y in L* e is 3a,2 + 0(e + 6) , so 
that all the coefficients of L* e have the same properties as the coefficients of Lq , where 
Lew = f/h is the equation obtained from L e u — f by setting -u = /wu. Therefore if 
g G L 2 (il), the problems 

L e w = g/h in fi, = 0, 

L* d w = g/h in fi, w x \ d n = 0, 

have weak solutions of the form w = v/h, where v G H^' 2 \Q) , i^i 1 ' 2 - 1 (SI) respectively. 
We then obtain 

Corollary 2.1. For eac/i g G L 2 {Q) there exists a weak solution v G H^ 1,2 ^ (f2) , 
#i 1,2) (fi) o/ (2.5), (2.6) respectively. 



3. Linear Regularity 



The purpose of this section is to establish the regularity, in X , of weak solutions 
to problem (2.3) for a particular choice of the right-hand side, /. This shall be 
accomplished by establishing the uniqueness of weak solutions to problems (2.3) and 
(2.4) in L 2 (£l), and then applying a boot-strap argument. 

In order to obtain the uniqueness of weak solutions, we will utilize the notion of 
a strong solution, in particular, for first order systems. The definition of a strong 
solution will be given below. We first introduce the notation and terminology that 
will be used for first order systems. Consider a boundary value problem 

SU = A 1 U X + A 2 U V + A 3 U = F in Q, U\ dn eN, (3.1) 
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where A 1 ,A 2 , A% are nx n matrices, U and F are n- vectors, and iV is a linear sub- 
space of the space of n-vector valued functions restricted to dfl. The corresponding 
adjoint problem is given by 



S*V = -A\V X - A* 2 V y + (A* 3 - d x A\ 



d y A*)V = G 



m 



n, V\ 9n eN\ 



where A* denotes the transpose of A i: and N* is the orthogonal complement of 
AN, where A is the matrix defined on dfl by A]U\ + A 2 n 2 , and (ni,n 2 ) is the unit 
outward normal to dfl. 

Let F G L 2 (Q). The notion of a weak solution to problem (3.1) is similar to the 
definition given in section §2 for single equations. That is, U G L 2 (Q) is said to be 
a weak solution of (3.1) whenever 

(S*V,U) = (V,F), 

for every V G C°°(fl) with bounded support and such that V\gn G N* . We now give 
the definition of a strong solution. 

Definition 3.1. U G L 2 (Q) is a strong solution of (3.1) if there exists a sequence 
{Uk]kLi, such that Uk G C°°(Vt) with bounded support, Uk\an £ N , and 



U k -U\\->0, || SUk — F ||— >■ 0, 



as 



k — > oo. 



Clearly, a strong solution is a weak solution. Moreover, using techniques developed by 
Friedrichs [2] and Lax-Phillips [14], G. Peyser [21] has obtained the following converse 
statement. 

Theorem 3.1 (Identity of Weak and Strong Solutions). Let the following 
conditions on the operator S and the boundary space N be satisfied: 

i) the matrix A is of constant rank in a neighborhood of the boundary, 

ii) N is of constant dimension at each point of the boundary, 
Hi) N contains the nullspace of A. 

Then a weak solution U G L 2 (Q) of (3.1) is also a strong solution. 

Note that for our particular domain, A = A]U\ , so that condition (i) is equivalent 
to Ai having constant rank in a neighborhood of dQ. 

With the aim of applying theorem 3.1, we shall transform problems (2.3), (2.4), 
(2.5), and (2.6) into the setting of first order systems. Let f,g G L 2 (Q) be the right- 
hand sides of (2.3), (2.4) and (2.5), (2.6) respectively, and define A\, A\, A 2 , A 2 , 
A 3 , A 3 , F, and G by 



A 1 = A 1 = 



( -9 a n \ 

1 





\ 1 / 



/ a 22 \ 




Ao = Ao = 






\ o o o 





1 





1 


0/ 
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Define boundary spaces N\ and N 2 by 

iVi = {(ui, . . .,u 5 )\ d n | u 5 \ dn = 0}, 

N 2 = {(u!, . . .,u 5 )\gn I {-Oui + a n u 3 )\ dn = 0}. 
Furthermore, define boundary value problems 

SqU = AiU x + A 2 U y + A 3 U = F in Q, U\ dn G N u (3.2) 

S e U = F in ft, U\ 9n e N 2 , (3.3) 

= + A 2 V y + A 3 V = G in Q, G N u (3.4) 

^\/ = G in ft, V\ dn eN 2 . (3.5) 

We now show that the weak solutions of (2.3), (2.4), (2.5), and (2.6) given by the- 
orem 2.1 and corollary 2.1 are also weak solutions of (3.2), (3.3), (3.4), and (3.5) 
respectively. 

Lemma 3.1. Let u G H^' 2 \Q) , H { }' 2) (9) be a weak solution of (2.3), (2.4) re- 
spectively, then U = (u xyy ,Uyy,u x ,Uy,u) G L? (ft) is a weak solution of (3.2), (3.3) 
respectively. Similarly, if v G H^ l,2 \fl) , H^' 2 \fi) is a weak solution of (2.5), (2.6) 
respectively, then V = {v xy y,v y y,v x ,v y ,v) G L 2 (ft) is a weak solution of (3.4), (3.5) 
respectively. 

Proof. Let u G H^ 2 \Q) be a weak solution of problem (2.3). We will show that 



/ / u*s;v = [ [ 

J Jn J Jn 



F*V 



(3.6) 



for all V G C°°(ft) with bounded support such that V\qq G N* , where 

N* = {(vi, ■ • - ,U5)|an I «i|an = ^lan = 0}. 
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A calculation shows that 

/ / U*SgV = (0u xyy - anu x )d x vi - a 22 u y d y v 1 - (ud x v 2 + u x v 2 ) 

J Jn J Jn 

+[(a 1 - d x a u )u x + (a 2 - d y a 22 )u y + au]v 1 (3.7) 

-(udyV 3 + UyV 3 ) - {UydyV 4 + UyyV^) ~ (UyyO^ + U X yyV^j. 

Since V\an G N* and u e H^' 2 \^l) is a weak solution of (2.3), we can integrate by 
parts to obtain 

f [ u*s;v = [ [ ul; Vi = I I f Vl = I [ f*v, 

J Jn J Jn J Jn J Jn 

showing that U is a weak solution of (3.2). 

Let u G H x ' (f2) be a weak solution of (2.4). We now show that (3.6) holds for 
all V G C°°(Q) with bounded support such that V\gsi G N£ , where 

^2 = {( v i, • • • ,^5) | an I v 2 \aci = v 5 \ dn = 0}. 
From (3.7) it follows that, 

/ / U*SgV = (0u xyy - a ll u x )d x v l - a 22 u y d y v x (3.8) 

J Jn J Jn 

+[(a 1 - d x a u )u x + (a 2 - d y a 22 )u y + au]vi. 

In order to integrate by parts we construct an approximating sequence for 
such that v\ e C™(Q) and 

II v i ~ v i II + II ®xVi ~ v i \\— > as k — > 00. 

Take a sequence {vk}^ =1 C C°°(f2) with the property that || Vk — d x v\ \\— > as 
k — > 00, and define 

/X 
v k (t,y)dt + v 1 (-x ,y). 
■xn 



Then since 



we have 



(wf-wi) 2 = (f d t (v k 1 (t,y)-v 1 (t,y))dt) 2 

J — xo 

< 2x I (d t v k 1 (t,y)-d t v 1 (t,y)) 2 dt, 

J —xo 

[ [(vt-vj 2 < 44 f [ (d x v k - d xVl ) 2 
J Jn J Jn 

= Ax 2 I I (v k - d xVl ) 2 -> 0, 
J Jn 
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so that Vi satisfies the desired properties. Therefore, recalling that a x | 
d x o>u\dn = by (ii) of lemma 2.2, and using the fact that a is a weak solution 
of (2.4), we can integrate by parts in (3.8) to obtain 




U*S;V = Mm I I uV B v\ 




= lim H= f Vl = F*V, 



n J Jn J Jn 



showing that U is a weak solution of (3.3). Similar arguments show that if 

is a weak solution of (2.5), (2.6) respectively, then 
V — (v xyy ,Vyy,v x ,Vy,v) G L 2 (Q) is a weak solution of (3.4), (3.5) respectively. □ 

Now that the weak solutions of the previous section have been translated into the 
setting of first order systems, theorem 3.1 is applicable. As a result, we obtain 

Proposition 3.1. The weak solutions of problems (2.3) and (2.4), given by theorem 
2.1, are unique in L 2 (Q). 

Proof. Let u G H^' 2 \il) be a weak solution of problem (2.3) with / = 0, then 

(L* e w,u) = for all w G C°°(Ji). (3.9) 

We will show that u = in L 2 (VL). 

Let v G H^' 2 \Q) be the weak solution of (2.5) with g = u. Then by lemma 3.1 
, v yy , v x ,Vy, v) is a, weak solution of (3.4). We now show that the conditions 
of theorem 3.1 are satisfied for problem (3.4). Condition (ii) is immediately satisfied, 
and since a* : < —9 in a neighborhood of dVt , condition (i) is satisfied with A = rkA\ 
having the constant rank of 3. Furthermore, the nullspace of A is given by 

{(v!, . . .,v 5 )\asi | (-Ovi + auv 3 )\dci = v 2 \an = v 5 \ d n = 0}, 

which is contained in Ni so that condition (Hi) is satisfied. Therefore, we can 
apply theorem 3.1 to obtain an approximating sequence {V^}^ for V, such that 
V k G C°°(VL) with bounded support, Vk\an G Ni, and 

\\V k -V\\^0, \\ S e V k - G as k ->■ oo. (3.10) 
From (3.10) it follows that 

II u fc - '"xw 0, || v 2 k -v yy 0, || vl -v x 0, 

II u * - u » IK °> ll w fe- w IK°> 

and 

|| (-9d x vl + a* u d x v 3 k + a* 22 d y v£ + a\vl + a\v\ + a*v\) - u \\-+ 0. 
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Hence, 

(u, u) = lim / / [-6d x v\ + a* u d x vl + a* 22 d y v\ + a\v\ + a* 2 v\ + a*v 5 k }u 

= lim / / (evl-a^vl^-a^vluy + lial-d^l^vl+ial-dya^v^+a^u 
J Jn 

= (8 v xyy - a* n v x )u x - a* 22 v y u y + [(a* - c^a^)^ + (a* 2 - d y a 22 )v y + a*v]u. 
J Jn 

Let {v n }^ =l C C°°(£l) be an approximating sequence for v in iJ( 1,2 )(fi). Then 
integrating by parts and using (3.9), we obtain 



lim / / (L* e v n )u = 0. 



Similar arguments hold for problem (2.4). □ 

Having established the uniqueness of weak solutions, we are now ready to apply a 
boot-strap procedure to obtain higher regularity for problem (2.3) in the x-direction. 

Theorem 3.2. Let u and f be as in problem (2.3). Let s < r — 4 and f G H s (£l) 
be such that d x f\gn = for a < s — 1. If s = e(s) is sufficiently small, then for all 
a < s, d^u G H^ l,2 \Vt) when a is even, and d^u G -ffi 1,2 -* (SI) when a is odd. 

Proof. The case s = is given by theorem 2.1. Consider the case s = 1. Let w = u x 
and formally differentiate the equation L e u = f with respect to x: 

L±w := -6w xxyy + a u w xx + a 2 2W yy + (a x + d x a u )w x + a 2 w y + (a + d x a{)w 
= fx- u yy d x a 22 - u y d x a 2 - ud x a := f 1 . 

Observe that since d x au, d x a\ = 0(e) and both vanish outside X, the operator Li 
has the same existence and uniqueness properties as Lg . Furthermore, by restricting 
Leu = f to the boundary of Q and using u\qq = ai\gn = 0, we obtain the following 
ODE 

(-0u xxyy + a u u xx )\ dn = 0, (3.11) 

for which the only solution in L 2 (dVi) is u xx \qq = 0. Therefore, in the regular case 
w = u x satisfies problem (2.4) with Lq and / replaced by L\ and f\. 

Let u G H^ 1,2 \VL) be the weak solution of problem (2.3). We now show that 
u x G L 2 (Q) is a weak solution of (2.4) with Lg and / replaced by L\ and f\ G L 2 (Q) ; 
we denote this problem by (2.4) 1 . Let v G C£°(f2), then 

(u x , L\v) = -(u, {L\v) x ) = -(u, L*(v x )) + (u, L*(v x ) - {L\v) x ) 

= -(/, v x ) + (u, -v yy d x a 2 2+ v y [d x a 2 - 2d xy a 22 ] + v[-d x a- d xyy a 22 + d xy a 2 }) 
= (fx, v) + {-u yy d x a 22 - u y d x a 2 - ud x a, v) = (fi,v). 
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Therefore weak solution of (2.4) 1 , and by the uniqueness result, proposition 

' — ' ( 1 2) 

3.1, u x must coincide with the solution in H x ' '(Q) given by theorem 2.1. Hence 

u x eH^ 2 \n). 

We now consider the case s — 2. Let u> = and formally differentiate the 
equation Liw^ = f\ with respect to x: 

L 2 w := -6w xxyy + a u w xx + a 2 2W yy 

+(a 1 + 2d x a 11 )w x + a 2 w y + (a + 2d x a x + d xx a u )w 
= d x fi - u xyy d x a 22 - u xy d x a 2 - u x (d x a + d xx a{) := f 2 . 

Again, since d x au, d xx au, d x a\ = 0(e) and all three vanish outside X , the operator 
L 2 has the same existence and uniqueness properties as Lg, provided that e is suf- 
ficiently small. Also, when u is regular u xx \gn = from (3.11). Thus in the regular 
case w = u xx satisfies (2.3) with L e and / replaced by L 2 and f 2 G L 2 {VL); we denote 
this problem by (2.3) 2 . 

Let u G H( 1,2 \Q) be the weak solution of (2.3), then we know that u x G H^ ,2 \Q) . 
We now show that u xx G L 2 (Q) is a weak solution of (2.3) 2 . Note that Lgu = f in 
L 2 (Vi) and let v G C°°(il), then a calculation produces 

(u xx , L* 2 v) = (u xxyy , -Qv xx ) + (u XX) (a u v) xx ) + {u yy , (a 22 v) xx ) + (u y , (a 2 v) xx ) 
+(u x , [(ax + 2d x a u )v) xx ) + (u, [(a + 28^ + d xx a u )v] xx ) 
= (L e u, v xx ) + (f 2 - f xx , v) = (/, v xx ) + (f 2 - f xx , v) = (f 2 , v). 

By the uniqueness of weak solutions for problem (2.3) 2 , u xx must coincide with the 
solution in H^ 2 \Q). Thus u xx G H^ 2 \tt) . 

To obtain the regularity of higher order derivatives, we observe that the above 
procedure applied to Lgu = f holds for L 2 u xx = f 2 , since for a > 1 

d^o-ii\dn = d x a 22 \gQ = d^ai\QQ, = d^a\gn = 0, 

so that f 2 \on = 0. Therefore u xxx G H x 1,2 \Q) and u xxxx G H^ 1 ' 2 ^^}). Furthermore, 
we can continue this process until / and the coefficients of Lg run out of derivatives, 
as long as e is chosen sufficiently small depending on the size of s. □ 

We now prove regularity in the ^/-direction for the weak solution of problem (2.3). 
The following standard lemma concerning difference quotients will be needed. 

Lemma 3.2. Let w G L 2 (Q) have bounded support, and define 

w h = —(w(x, y + h) — w(x, y)). 

I i 

U || wh || < M where M is independent of h, then w G if( 0,1 )(r) for any compact 
fcO. Furthermore, if w G H^'^^Q) then || w h \\< M \\ w y \\ . 
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Theorem 3.3. Let the hypotheses of theorem 3.2 hold, then u G H S (X) . 

Proof. From theorem 3.2 we know that d x u G H^' 2 ^(Q) for < a < s. Therefore 
the following equality holds in L 2 (Q), 

Luyy := -9u xxyy + a 22 u yy = / - a u u xx - aiu x - a 2 u y - au := /. (3.12) 

Since \a 2 \ = 0(\y\) as \y\ — > oo, we do not necessarily know that / G H^ 0,1 \il) ; 
however, we do have / G H^ 0,1 \T) for any compact r C O. Fix a constant k > y 
and set w = v k u yy , where v k is given by (2.7). Then 

Lw h = {vjf - v k {y + h)u yy (x, y + h)a h 22 . (3.13) 

Since u G H^ 2 \Q), by multiplying (3.13) on both sides by w h and integrating by 
parts, we obtain 

|| ^ll + ll w h x \\<M 1 (\\ [y k f) h || +1), 
for some M 1 independent of h. By lemma 3.2 

II w h ll + H w h x ||<M 2 (||i/ fc / 11(0,1) +1), 

independent of h. Therefore w y , w xy G L 2 (X), which implies that dyU, 
d x dyU G L 2 (X). Furthermore, by differentiating Lqu = f with respect to x, 
a = 1, . . . , s — 3 times, the same procedure yields d^dyU G L 2 (X) . 

Proceeding by induction on /, assume that d^d^u G L 2 (X), a < s — /3, (5 < I, 
and 3 < I < s. Differentiate (3.12) with respect to y, I — 2 times: 

1-3 

LdyU = d l y 2 f - J2 dl(dya 22 d l y^uyy). (3.14) 

i=0 

Note that this equation holds in L 2 (Vt), and that the right-hand side is in if( 0,1 )(r) 
for any compact rcfi. Applying the method above yields d l y +1 u, d x d l y +l u G L 2 (X). 
Moreover, differentiating (3.14) with respect to x, a — 1, . . . , s — (I + 1) times, and 
applying the same procedure, yields d x d l y +l u G L 2 (X). The desired conclusion now 
follows by induction. □ 

4. The Moser Estimate 

Having established the existence of regular solutions to a small perturbation of the 
linearized equation for (1.5), we intend to apply a Nash-Moser type iteration proce- 
dure in the following section, to obtain a smooth solution of (1.5) in a subdomain of 
X which contains the origin. In the current section, we shall make preparations for 
the Nash-Moser procedure by establishing a certain a priori estimate. This estimate, 
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referred to as the Moser estimate, will establish the dependence of the solution u 
of (2.3), on the coefficients of L e as well as on the right-hand side, /. The Moser 
estimate that we seek has the form 

II u \\h'(X)< C,(\\ f \\h'(x) +A s+S() II / \\h*{x)), (4-1) 

where 

A-s+so = II a ij \\h s + s o(X) + II a i \\h s+s o(X) + II a \\h s + s o(X) 

for some s > , and C s is a constant independent of e and 6 . 

Estimate (4.1) will first be established in the coordinates {C,,f]), which we have 
been denoting by (x, y) for convenience, and later converted into the original coordi- 
nates (x, y) of the introduction. We will need the following preliminary lemmas. The 
first is a modification of lemma 2.3, and the second contains standard consequences 
of the interpolation inequalities for Sobolev spaces. 

Lemma 4.1. Let w G H^^Q) (or H X 2 ' 2 \Q)) be such that yw G L 2 (Q) , and let 
Pi = e0pi, P2 = sp2, P3 = eps, where pi G Cf{X), i = 1,2,3. Then for e and 9 
sufficiently small, there exists a constant M , independent of e and 6 , such that 

II HI + II w y \\< M || piw xyy + p 2 w x + p 3 w + Lqw II . 

Proof. Assume temporarily that w G C°°(f2) (or C£>(Q)). The properties of p 2 and 
p 3 guarantee that lemma 2.3 holds for the operator p 2 d x + p 3 + L 6 . Therefore 

(Aw + Cw y + Dwyy, p 2 w x + p^w + L qw) > (4.2) 

Cl[|| W f + || Wy f +0(\\ W X || 2 + || W X y || 2 + || Wyy f)], 

where A,C,D, and C\ were given in lemma 2.3. Furthermore, integrating by parts 
yields 

(AW + CWy + DWyy,PlW X yy) = ^ ^ [ ~ ^ (Dp l ) J W^ y + [~ ] W X y Wyy 

+ [\(Cpi) xy + l -(Ap l ) x }w 2 y (4.3) 

+ [(Ap l )y]W X Wy + [-^(Ap^^yjW 2 

All the boundary integrals vanish since p\ G C^°(X). Moreover, the properties of 
Pi guarantee that by choosing e and 6 sufficiently small, we obtain the following by 
adding (4.2) and (4.3), 

(AW + CWy + DWyy, P\W X yy + P 2 W x + J9 3 W + LqW) > 

2 i II „.. 1 1 2 i q/\\ „.. 1 1 2 i || „.. 1 1 2 i || „.. ||2M 



Cl[|| W || 2 + || Wy || 2 +^(|| W X || 2 + || W X y || 2 + || W 



yy 
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Then an application of Cauchy's inequality, and the use of an approximating se- 
quence {wk}^ =1 , as was constructed in proposition 2.1, removes the assumption that 
w G C°°(Q) (or and completes the proof. □ 

Lemma 4.2 [24]. Let u,v e H S (X). 

i) IfQ<i<j<s, then there exists a constant A4ij tS such that 



3—3 J— I 



II u \\w{X)< Mij, 8 || u \\ s H i\ x) \\ u \\ s H s\ x) ■ 

ii) If a and (3 are multi-indices such that \a\ + \(3\ = s, then there exists a 
constant A4 S such that 

|| d a ud f) v || L 2(x)< M s {\u\ L ™( X ) || v \\h>(x) + II « \\h*(x) \v\l°°{x))- 

in) Let T C M N be compact and contain the origin, and let G G C°°(T). If 
u G H S+2 (X, T) and \\ u \\h 2 (x)< C for some fixed C, then there exists a constant 
M. s such that 

|| Gou \\h*(x)< Vol(X)\G(p)\+Ms II u \\ H '+*{x) ■ 

Estimate (4.1) will be established by induction on s, and we begin by estimating 
the a; -derivatives. Let || • || s ,x denote || • and | • |oo denote | • \l°°(x)- 

Proposition 4.1. Let u and f be as in theorem 3.2. If e = e(s) is sufficiently small, 
then 

|| d s x u || + || d s x u y \\< C s (\\ f \\ a + || u \\ s -i,x +A s +2 || / Ikx), 
for s < r — 6 , where C s is independent of e and 6 , and 

A s +2 = || 0>ij ||s+2,X + || a i ||s+2,X + || a ||s+2,X • 

Proof. We proceed by induction on s. The case s = follows from lemma 2.3. 
Differentiate Lgu = f s-times with respect to x and put w = d x u, then 

- 9w xxyy + a n w xx + a 2 2W yy + (ai + sd x au)w x + a 2 w y + a s w (4.4) 

s— 1 s— 1 S — 1 

i=0 i=0 i=0 

:=/. 

where a s = a + sd x a\ + s ^~^ ^an . A calculation shows that 

S_1 • sis — 1) 
^2d x (d x a 22 d s x ~ 1 ~ t Uy y ) = sd x a 22 d s ~ l Uyy^ d 2 x a 22 d x ~ 2 u yy 
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+ EE(*)^ +I ^r 1 --' 



Uyy. 



i=2 j=2 



Note that the term d s x ~ x u yy contains too many derivatives. However, since 
a 2 2 = 1 + 0(e), we can solve for d s x r l u yy in (4.4) with s replaced by s — 1 to 
obtain a more manageable expression: 

dx~ lu yy = —[Qwxyy - a n w x - (ai + sd x a n )w - a 2 d s x ~~ 1 u y - a s ^ x d s x r l u + / s _i]. 

022 

Substituting back into (4.4), we have 

s9d x a 22 , / a sa u d x a 22 . , , sd x a 22 , 

+ (sd x an )w x + (a s - a (ai - sd x an))iu + L e u? 

^22 a 22 a 22 

= 3/ - S -^dla 22 d^u yy - ( ) )di +1 a 22 d:-^u y y-^dUd x a 2 d s x -^Uy) 

i=2 j=2 ^ ' i=0 

- V di&as^d^u) + S -^[a 2 d s - l u y + a^diT'u - /U] 
i=o ° 22 

If £ = e(s) and are sufficiently small, we can apply lemma 4.1 to obtain 

\\d s x u\\ + \\d s x u y \\<M\\J s \\. (4.5) 

We now estimate each term of f s . Using lemma 4.2 (ii) , lemma 2.2 (Hi) , and the 
fact that d x a 22 vanishes outside of X , produces 

ii EE ( !•) ^^r 1 -^ ii=ii EE ( •) ^^r 1 -^ ik* 

i=2 j=2 \ / j =2 j= 2 W / 

< Mi(|^a 22 |oo || « || s -i,x + || dla 22 \\ s -i,x Hoc) 

< M[(\\ U \\s-1,X + || a 22 \\s+2,X \u\oo). 

A calculation shows that 

s— 1 s — 1 i / • \ 

£ dl(d x a 2 d^u y ) = 8d x a 2 d'-\ + E E ( ' 
i=o i=i j=i V ^ / 

Then using the same procedure as above, we have 

s-l 

|| ^^(c^c^X) ll< M 2 || ^r 1 ^ II +M' 2 (\\ u \\ s ^ x + || a 2 || s+2 ,x |«U). 

i=0 
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Furthermore, the following estimates are obtained in the same way: 

|| f^di&a.-^- 1 -**) || + || ^f^^flS-^) || 
< M 3 (|| it || s _i ;X +(|| a \\ 3+ 2,x + II «i || s +2,x + || an || s +2,x)Moo), 

and 

sd a s ~ 2 

|| - ^2d l x (d x a 2 d x ~ 2 ~ l u y ) ||< M 4 (|| u || s _ 1>x + || a 2 || s+2 ,x |u|<x>)- 

«22 J; 
1=0 

Also, since 

s— 2 s— 2 i / • \ 

£ ^(9xa22^- 2 - i « w ) = (s - l)d x a 22 d s x - 2 u yy + ££(!•) 2 

i=0 i=l j=l ^ ' 

and c^g^ = 0(e), we find that 

|| ^^(d^d^uj \\< es 2 M 5 \\ d^u yy || ,* 

a 22 n 
i=0 

+ M' 5 (\\ u \\ s -i,x + || a 22 \\s+2,X Hoc), 

where M 5 is independent of e and s. 

Summing the above estimates produces: 

II Is ||<M 6 (|| / || s + || u || s -i,x+|| d s -\ \\+es 2 || d s ~ 2 u yy \\ ,x +A a+2 |u|oo)- (4.6) 

Therefore, if we estimate || d x ~ 2 u yy || ,x appropriately and show that 

Moo < M 7 || / \\ 2tX , 

the proof will be complete by induction. 

We now estimate || d x ~ 2 u yy \\o,x ■ Differentiate the equation, 

Lu yy := -6u xxyy + a 22 u yy = f - a u u xx - a x u x - a 2 u y - au := g, 
with respect to x (s — 2) -times, then 

s-3 

Ld s x ~ 2 u yy = d s x - 2 g - dKd^d^Uyy) := g s ^ 2 . 

i=0 

Multiply the above equation by d x ~ 2 u yy and integrate by parts in X to obtain, 

II 9 S x~ 2u yy \\o,x< M 8 \\ g s -2 \\o,x ■ 



33 



We now estimate || g s _ 2 \\o,x- Using the same methods as above, we have 

s-3 

|| d s x ~ 2 (a 1 u x + a 2 u y + au) + ^ d x {d x a 22 d s x ^~ l Uy y ) \\o,x< M 9 (\\ u \\ s -i,x +J\ s +2\u\oc)- 
Furthermore, 

s-2 , x 

d x - 2 (a n u xx ) = a n d s x u + ^ ( S ~ j ^11^^; 

i=i ^ ' 



thus, 

|| d x ~ 2 (a u u xx ) \\o,x< Afio(|| d> llo,x + || u \\ s -\,x +A s+2 
It follows that 



u 



II d s x 2 u yy || 0l x< M n (|| d> || 0l x + II « +A s + 2 |m|oo). (4.7) 

The coefficient of || d x ~ 2 u yy \\o,x in (4.6) is es 2 Mq. If £ = e(s) is chosen sufficiently 
small so that es 2 MM 6 M n < |, we can then bring es 2 MM 6 Mn || <9*w || ,x from (4.7) 
to the left-hand side of (4.5), so that by induction on s 

II d s x u || + || d s x u y \\< M' e (\\ f \\ s + || u || S _ 1)X +A s+2 (|m|oo+ II / || 2 ,x)). 

We now estimate |u|oo to complete the proof. The above methods can be used to 
show that 

II u h,x< M 12 || / || 2 ,x • 

Then by the Sobolev lemma, 

i«ioo < m 13 11 u i| 2iX < m; 3 11 / || 2>X . □ 

We now estimate the remaining derivatives. 
Proposition 4.2. Let u, f , s, and e be as in proposition 4.1. Then 

II d^u || ,x< C s (|| / \\ s , x + || u || S _ 1)X +A s+2 || / || 2iX ), 

for a + (3 < s, where C s is independent of e and 6 . 

Proof. The cases f3 = 0,1,2 follow from (4.7) and proposition 4.1. We proceed by 
induction on (3. Assume that the desired estimate holds for < a < s — (3, and 
< (3 < k — 1 , for some k < s. 
Differentiate the equation, 

LUyy := -6U XX yy + a 22 Uyy = f ~ CLuUxx ~ 0\U X - 2 Uy - ttU 
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with respect to d^d y 2 where < a < s — k, then 

a-l fc-3 

Ld a x d k u = d a x d k y - 2 g - d^dHd^dr'-'uyy) - dl(d y a 22 d k y - 3 - l d:u yy ) 

i=0 i=0 

'■= 9a,k-2- 

Multiply the above equation by d x d k u, and integrate by parts in X to obtain 

II & y u \\ ,x< M || 9a,k-2 \\0,X ■ 
We now estimate || g a ,k-2 \\o,x ■ Using lemma 4.2 (ii), we have 

|| d:d k y ~ 2 (a n u xx ) \\ ,x < M^ll d: +2 d k -\ \\o,x + £ H d^ y andr p d k ^u xx \\ ,x) 

p<a, q<fc — 2 
(p,g)^(0,0) 

< Mi(|| 9 x +2 d k ~ 2 u \\o,x +|an| c i ( x) IM| s -i,x + ||an IUxMoo) 

< Mi (|| d a x +2 d k ~ 2 u || ,x + || u \\ s -i,x +K+2 || / || 2 ,x). 

Furthermore, if a < s — k then || d x +2 d k ~ 2 u ||o,x<|| u || s -i,x, and if a = s — k the 
induction assumption implies that 

II d: +2 d k ~ 2 U || ,X< M 2 (|| / \\ 8tX + || U \\ s _ 1>x +^s + 2 || / h,x). 

Thus, 

|| d:d k y - 2 (a n u xx ) || ,x< M 3 (|| / || s ,x + || u h-^x +A s+2 || / || 2 ,x). 

Moreover, the methods of proposition 4.1 may be used to estimate the remaining 
terms of || g a ^2 \\o,x by 

M 4 (|| U \\ s ^x +^s + 2 II / || 2 ,x). 

The desired conclusion now follows by combining the above estimates. □ 

From proposition 4.2, we obtain the following Moser estimate by induction on s. 

Theorem 4.1. Let u and f be as in theorem 3.2. If e = e(s) is sufficiently small, 
then 

II u \\ SjX < C s (\\ f \\ S:X +A s +2 || / h,x), 

for s < r — 6 , where C s is independent of e and 9 . 

The estimate of theorem 4.1 is in terms of the variables (£,77) of lemma 2.2, 
which we have been denoting by (x, y) for convenience. We now swap notation and 
denote the original variables of (2.1) by (x,y), and the change of variables by (£,77). 
Furthermore, let || ■ || s and || • || a denote the H S (X) norm with respect to the 
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variables (x,y) and (£,77) respectively. Similarly for A s and A^. We now obtain 
the analogue of theorem 4.1 with respect to the variables (x,y). We will need the 
following lemma. 

Lemma 4.3. If s — e(s) is sufficiently small, then 

|| 6r ||s< Cs (|| a 12 \\s+3 + || a 22 \\s+b), 

for s < r — 7 , where C s is independent of e and 6 . 

Proof. We prove the estimate by induction on s. The case s = follows from the 
estimate, 

0<M 1 < 16,1 < M 2 , 

obtained in the proof of lemma 2.2. Now assume that the estimate holds for s — 1. 
We first estimate the ^-derivatives. Differentiate the equation 

(—)(&)* + (&), = -(—)*& (4-8) 

a 22 a 22 

with respect to x s-times to obtain, 

(—)(%). + (%)„ = -d s x [(— u»i - yr i u := h s . 

0-22 0-22 n &22 

Then estimating d s J^ x along the characteristics of (4.8) as in the proof of lemma 2.2, 
we have 

\d S xix\c*(x) ^ M 3\h s \ C o(x)- 

Recalling that 012 = 0(e) , and using the analogue of lemma 4.2 (ii) for C S (X) norms 
in the same way that the Sobolev version was used in proposition 4.1, produces 

\h s \ c °(x) < £(s + l)^4|%| c o(x) 

+M 4 (\( — )xx\cO(X)\^\c s - 1 (X) + \ {~)xx\c^(X)\ix\c»(X))- 
022 0,22 

Therefore, if e is small enough to guarantee that e(s + 1)M^M4 < |, we can bring 
e(s + 1) M 3 M 4 \d^ x \ c o^ to the left-hand side: 

\dx£x\c°(x) < M 5(\Cx\c- 1 (x) + I" - lc*+i(x))- (4-9) 

a 22 

We now estimate the remaining derivatives. Assume that 

\%tft*\co<X) < Med&c-w + |^| cs+1( x)) (4.10) 

°22 
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for all < a < s - /3, < (3 < s - 1 . The case (3 = is given by (4.9). Differentiate 
(4.8) with respect to <9" -1 <9^ to obtain 

sr 1 ^. = )(^x)x] uj - ^E^cc^)^- 1 -^). 

1=0 

Using assumption (4.10) on the first term on the right-hand side, and applying lemma 
4.2 (ii) to the remaining terms, we find 

\&x 1 dy +1 £x\co(X) < ^{\ix\c^(X) + \~\c a + 1 (X))- 

0,22 

Thus, by induction on (3, estimate (4.10) holds for all < a < s — (3, < (3 < s . 
By induction on s, (4.10) implies that 

I6z|c s (x) — M s \ — \ C s+i(x)- 



Then the Sobolev lemma gives 



^22 



6r ||s< Mg || — || s+3 
a 22 



Moreover, by lemma 4.2 (ii) and (Hi) we have 

|| ||s+3 < A^lo( |^12 |oo || ||s+3 + || a 12 \\s+3 | |oo) 

022 0-22 O22 

< Mn (|| a 22 || s +5 + || 012 || s+s)- □ 

Theorem 4.2. Let u and f be as in theorem 3.2. If e = e(s) is sufficiently small, 
then 

II u \\ s < C s (\\ f \\ s +A s+ n || / || 2 ), 
for s < r — 13, where C s is independent of e and 9 . 
Proof. Let a be a multi-index with \a\ < s. A calculation shows that 

||^> || < Afx || G i d ln u I 

h\<s 

where G 7 are polynomials in the variables x^ 1 = £ x , d^x^, and dj 2 x v , such that 
J2i < s — I7I for each term of G 7 . Then using lemma 4.2 (ii) and (Hi), we find 
that 

II d a £ , y u || < M 2 (|| u ||1 +(|| Xi : \\' 8+2 + || x v ||l +2 )K|oo). 

Similarly, 

II dl ri u || < M 3 (|| u \\ s +(|| & || s +2 + || iy || s+2 )k|oo)- (4.11) 
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Then by theorem 4.1 and the Sobolev lemma, we have 

II ^ y u || < M 4 (|| / Hi +Al +2 || / ||' 2 ) + M 4 (|| Xi ||l +2 + || x v ||l +2 ) || / || 2 . (4.12) 

We now estimate the terms on the right-hand side of (4.12). Use lemma 4.2 (ii) , 
(in), and (4.11) to obtain 

M M 7 m^ii' I. . 1 1 / 

II X Z lls+2 = II T~ lls+2 ^ M 5 || £ x || s+4 

< M 6 (|| £x ||s+4 +(|| 6r l|s+6 + || £,y ||s+6) |<vr|oo) 

< M 7 (\\ £ x \\ s+6 +\\ — £ x || s+6 ) 

< M 8 (|| Oi2 || s +9 + || 022 lU+n) - 

Similarly, 

II X V lls+2 = II 7~~ L+2 < ^9(11 fl 12 ||s+7 + || ^22 || s+s) - 
Si 

Furthermore, 

H/li; < M 10 (||/|| S +(|U, || S +2+|UJ| S +2)|/|oo) 

< M 11 (\\ f \\ s +(\\ a 12 || 

s+5 + || a 22 ||s+7j || / Ih), 



and hence 
Also, 



2 < M 12 (|| a 12 || 7 + || a 22 || 9 ) || / || 2 < M 13 



0>ij || s+ 2 < ^14(|| % ||s+2 +(|| £x ||s+4 + || £2/ ||s+4)| a ij|oo) 
< -^15 (|| °ij ||s+2 + || a 12 \\s+7 + || a 22 ||s+9,)> 



so that 

Al +2 < M 16 A s+9 . 

Therefore, using the above estimates and summing over all |er| < s, (4.12) produces 

II u \\ s < M 17 (\\ f \\ s +A S+11 || / || 2 ). □ 

5. The Nash-Moser Procedure 

In this section we will modify the Nash-Moser iteration procedure to obtain a solution 
of 

$H = in Xoo, (5.1) 
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where C X is a sufficiently small neighborhood of the origin that will be defined 
below. In order to accommodate the requirement (theorem 3.2) that d%f\dn = 0, 
a < s — 1 , we will cut off the right-hand side of the modified linearized equation, 

L e u = f, 

near dX at each iteration, and then estimate the error in a smaller domain at the next 
step. Furthermore, the constant 9 will be chosen sufficiently small at each iteration, 
to guarantee that the procedure converges. 

Let /i > 5. Define a sequence of domains {X n }^ =1 by 

n-1 

X 1 = X, X n = (l-^/i- l )X, 

where AX = {Xx \ x <E X}. Then X^ = (1 - ~z^)X . In addition, let | < r < 2 and 

define /j, n = fi rn+n ° , where hq > will be chosen sufficiently large. 

We now construct smoothing operators on L 2 (X n ). Fix ip G C^°(M. 2 ) such that 
•0 = 1 in some neighborhood of the origin. Let tp(x) = J J R2 tjj(r])e 2mvmx dri be 
the inverse Fourier transform of ip ■ Then ip is a Schwartz function and satisfies 
/ J R 2 ip{x)dx = 1 , and J j R2 x a ip(x)dx = for any multi-index a ^ 0. If g G L 2 (R 2 ) 
and 7 > 1, we define the smoothing operators S' : L 2 (IR 2 ) — > H QO (R 2 ) by 

(S'g)(x) = -f 2 / ip(-f(x -y))g{y)dy. 
J Jr 2 

Then we have (see [22]), 

Lemma 5.1. Let a, b G Z> and g G H a (M. 2 ) , then 
i) II S^g ||ff&(]R2)< C a;6 || # ||j?a(R2), & < a, 

II s '~,9 \\h»(ri)< C a ^ b - a || 5- ||jy«CR2), a<b, 
Hi) || 9 - S' 7 g || H 6(R2)< C ai6 7 6 " a || 5 ||h«(r2), 6 < a. 

To complete the construction, we also need the following extension theorem. 

Theorem 5.1 [23]. Let D be a bounded convex domain in IR 2 with Lipschitz smooth 
boundary. Then there exists a linear operator Td : L 2 (D) — > L 2 (IR 2 ) such that: 

i) T D (g)\ D = g, 

ii) T D : H a (D) — > H a (M. 2 ) continuously for each a G Z> . 

To obtain smoothing operators on X n , S n : L 2 (X n ) — > H°°(X n ), we set 
'S'nfi' = (»S' Mn 7x n fi , )U„ • Furthermore, it is clear that the corresponding results of lemma 
5.1 hold for each S n . 

We now set up the iteration procedure. A sequence of functions {w n }™ =1 will 
be shown to converge to a solution of (5.1), and shall be defined inductively as fol- 
lows. Set w± — and suppose that Wj, j < n, are already defined in Xj, then set 
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w n +i = w n + S n u n in X n+ i , where u n is defined in X n and will be specified below. 
Set /„ = — $(u>„) in X n , and let n be a C°° cut off function given by 



l if x e 

if x G X - X n , 



such that 
Let 



\<t>n\c°{X n ) < M sf i Sn . 

L ( w n) = ^2 a ij( w n)dij + ^2 ai(w n )di + a{w n ) 



1,3 i 

denote the linearization of evaluated at w n , and let {9 n }™ =1 be a sequence of 

positive numbers tending towards zero that will be specified later. Then define u n in 
X n by u n = v n \x n , where v n is the solution of 

Le n (Wn)V n = <Pnfn "1 X, 

given by theorem 2.1. Since fi > 5 we have |X C X^. Therefore, it follows from 
the definition of $(u>) in (1.5) that the coefficients of L 8n (w n ) are well-defined in all 
of X , even though w n is only defined in X n . 

For simplicity we denote the Sobolev norms || • ||/fs(x n ) by || • ||™, and the C s (X n ) 
norms by | • |™. Let s* G Z> be fixed such that $(0) G H S *(X), and define 

a = n(n+l)T- { - n+l+na \ 5 



r-1 



The convergence of the sequence {w n }™ =1 to a solution of (5.1) will follow from the 
following four statements. Each will be proven by induction on j , for some constants 
C\ , C2, and C3 independent of j and dependent on ji and s* . We shall require that 
s < s* - 18 - 25 - £p and s* > 22 + 25 + ^jz_ . 

T ■ II ?/>• IP < n us+s II f-, II 1 

i J • II W J lls+15^ H'j II Jl lis*— 15 



-15 



u i- II lis < ^1^-1 II Ji L, 

III, : || £ ||j< C ^ s - s * +1S+25) HMIU5 

IV, : || «;,- ||{ 4 <C 3 

To start the induction process observe that Ii, II 1, and IV 1 are trivial, and that 
IIIi holds if we set C 2 = H\. Now assume that I,, . . .,1V, hold for 1 < j < n. The 
next four propositions will prove the induction step. Note that the coefficients of 
L(wj) satisfy the conditions placed on (2.1) with r = — 2. Therefore, the results 
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of the previous sections apply to Lg.(wj), 1 < j < n, as long as e(s*) and 9j are 
sufficiently small and s < s* — 15. 

Proposition 5.1. // s < s* — 15 and /i(s*) is sufficiently large, then 

|| w n+ i || s+15 

< ACi II A II s» — 15 ' 

Proof. We have 

II lls+15 — II W « IL+15 + II 'S'n^ri || s+15 • 

Furthermore, by theorem 4.2 and lemma 4.2 (Hi), 



S n u n \\: +15 < M^\\u 



n Us 



< M 2 ^ 5 (||0 n / n ||:+||«; n ||: +15 ||0 n / n ||^). 
Using lemma 4.2 (ii), we obtain 

HnfnW: < M 3 (ii/ n |i: + ii0 n |i:i/„p 

< M 4 (||/ B ||»+||^ B ||»||/ B ||5) 

< M 5/ /" || / n II? . 

Moreover, by definition of /„ and lemma 4.2 (Hi) 

|| fn \\s< M e(\\ fl ||?„-15 + II W n ||"+4)> (5-2) 

so that, 

|| 0n/n C< M 7/ /"(|| /l||I_ 15 +|| W n C +4 ). 

Similarly, using IV „ 

II 0n/„ ||J< M 7/ , 2 "(|| fi ||:._ 16 + || «;„ ||J) < M 8/ i 2 ". 

We now have 

|| S„U n ||: +15 < M 9 ^V Sn (ll fl llL-15 + II ^ c +15 ). 

Therefore, 

|| Wn+1 ||H 6 < 2M 9 ^V"(II fl IIL-15 + II ||: +15 ) 

< ^v 2sn (ii /i iiL-15 + ii «>» ur+j. 

where the last inequality holds if /i is chosen so large that 2M 9/ u~ 1 < 1. It follows 
that 



Wn + 1 \\%h< (tl^ M )M 10 HAUL 



Is*— 15' 

i=l 
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where 

n-l 

m w = i + ^r 1 v 2s + + n V 2 " < 2, 

i=i 

if ji is large. Hence 

II «Wi ||S 6 < 2^^ 1 )^^ +1+ "°- 1+ "°) || /1 ||^ 15 
<r II f II 1 

— Pn+l II il lls»-15' 

where a = n(n + i) r -(«+ 1 +«o) and 5 = 

Proposition 5.2. // s < — 20 — 25 and n (s*) is sufficiently large, then 

II n, \\ n < O ,,">- 1 (s-s*+18+2(5) 11 r ||1 

|| «n L< || /1 || St _i 5 , 

where C\ depends on fi and s* . 
Proof. By theorem 4.2 

|| «n HI _ 15 < Mid) 0n/n ||™ _ 15 + || W n <f) n f n ||»), 

where Mi depends only on s*. By lemma 4.2 (ii), (5.2), and I n 

II <Pnfn ||?._15 < M 2 (|| / n ||« „ 15 + || n ||^_ 15 || f n ||») 

< m 3 (i + ^- 15 )"K^- 26 ) +5 ||/ 1 || s \_ 15 

< M 4/ i^ +5 11 M|^ 15 , 
where M3 depends only on s*. Similarly, III n yields 

II 0n/„ ||5< M 5 C 2/ ir T " (20 ^* +25) II fl IIL-15 • 

Therefore, for some constant M 6 depending on \i and , we have 



11-15 < M 6 (^^+<(-- i5 Hvr T " 1(20 ^ +25) )iiAiii-i 

< 2M 6 ^-+ 5 || /x ||i 



since s* > 20 + 25. Furthermore, lemma 2.3 and III n produce 

II u n ||»< M 7 || /„ ||»< M 7 C 2/ .r 1(18 -" +25) II fi Ill-is • 
Then applying lemma 4.2 (i), we find 

IK II: < M 8 (\\ Un ii^-^q^iii^)^ 

<T A/f , -- 1 (18-^+2^)(l-^)+(2s^+5)(^ T 3) I, .,! 

S M 9 /i„ || /1 || s ,_ 15 

< M 9 ^- 1 C-+ 18 + M ) || /1 IIL-15, 
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if a is sufficiently small. Note that a may be made arbitrarily small by choosing n 
sufficiently large. We then set C 2 = M 9 to obtain the desired result. □ 

Proposition 5.3. // s < s* - 18 - 25 - ^ , s* > 22 + 25 + ^_ } no ( s ^) a nd fi(s*) 
are sufficiently large, and e(s*) sufficiently small, then 
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II ^ n ,. T 1 (s-s*+18+2<5) 11 j. Mi 

|| Jn+1 || s S <^2/Vf-l || Jl || s »- 

Proof. Expanding $(u> n+ i) in a Taylor series yields, 

/n+1 — fn ~ L{w n )S n U n + Q n = f n — O n {S n U n ) m ^ — Lg n (w n )S n U n + Q n , 

where (£,77) are the change of variables given in section §2 by 

ai 2 (wn)tx + a 2 2(w n )i y = in X, £(x, 0) = x, f (=bc , y) = ±x , rj = y, 
and where Q n is the quadratic error term given by 

Qn= [ (t-l)d*®(w n + tS n U n )dt. 

Jo 

Since L 6n (w n )u n = f n in X n+1 , we have 

fn+i = Le n {w n ){u n — S n u n ) — 9 n (S n u n ) m ^ + Q n , (5.4) 

in X n+1 . 

Each term of (5.4) shall be estimated separately. First note that 9 n may be chosen 
sufficiently small to guarantee that, 

II (Q ,, \ \\ n + l <- ^ r i ,t- 1 (s-s,+18+2S) H , Mi 

|| ^n\O n U n ) m ^ || s S 2^2^ n+ i || Jl L.-15 • 

We now estimate L dn (w n )(u n — S n u n ) . By lemma 4.2 and IV n , 

|| Le n (w n )(u n - S n u n ) \\ n s +1 < || L en (w n ){u n - S n u n ) \\ n s 

< Mx(|| u n - S n u n \\ n s+2 + || w n \\ n s+A \u n - S n u n \%) 
+0(8 n ) 

< M 2 (|| u n - S n u n \\ n s+2 + || w n \\ n s+i \\ u n - S n u n 
+0(6 n ) 

< M 3 « +17 -- || u n ||« _ 15 +pt 1 J- s * || w n \\: +4 \\ u n \\l _ 15 ) 
+O(0 n ). 



Furthermore, by (5.3) 



«n ||?.-16< Mitf*** II h IIL- 



15 
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If 6 n and a are sufficiently small and /x is sufficiently large, it follows that 

II T {on \(m Q 11 \ 11™+! <? A/f ., — l7..3s»<T+s— s*+17+5 , 3s»cr-s» + 17+2<5\ n /• nl 

II L e n {w„){u n - b n u n ) \\ s < M 5 /i +/i n J II Ji IU 



s* — 15 



/ l^y , t- 1 (s-s„ + 18+25) I, , nl 

S 2^2A t „+i || Jl || S ,-15 • 



We now estimate Q n . Apply lemma 4.2 (ii) to obtain, 

ii o ir +1 < ii o ii" 

|| Us _ || lis 

< f \\d P Hw n + tS n u n )d a (S n u n )d^S n u n )\\: dt 

J ° l«l,l/3|,H<2 

< [ J2 M 6 (\d^(w n + tS n u n )\Z\\d a (S n u n )dP(S n u n ) \\ n s 

JO uiall-l/n 



M,I/3|,W<2 

Then the Sobolev lemma and the interpolation inequality || w 2 || L 2< C \\ u llj^-i , show 
that 

II Qn WT 1 < f Yl M 7(H + ^n) II 2 (II SUn ||? +3 ) 2 

^ |p|<2 

+ ii^K+*s i n « n ) n:(ii^„ ii^) 2 )^. 

Furthermore, by lemma 4.2 (Hi), I„, IV„, and proposition 5.2, 

II QnC +1 < M 8 [(|| «; B ||?+^ || u n \\ n M || ^C) 2 + (|| «; B ||J +4 +^ || unl^dl u n \\lf] 

< (M B II A ||s,-i 5 )[(l + ^ r " ( ^* +22+25) )^ +2r " (s -" +18+25) 

+ (< (s - ii)+i +^ 1(s ^ +i8+M vr 1( - s * +22+25) ] ii /i iiL_ 15 

< (Mo || fl Hs.-lsVn S * + 18+2<5 || /i ||s.— 15 5 

since s < s* — 1 8 — 25 — and > 22 + 25 + . If e(s*) is sufficiently small to 
guarantee that M 10 || / x ||s <i _ 1 5< |C 2 , then 

ii n \\n+l^ 1 r> , ,t- 1 (s-s*+18+2S) || nl 

|| Vrt ||s ^ 2 < - / 2A t n+l || Jl Ms* — 15 ■ 

By combining the estimates for each term of (5.4) we obtain the desired result. □ 
Proposition 5.4. If n (s*) is sufficiently large, then 

|| UJn+l Hlz^ 1 ^ C3, 

where C3 depends on ji and s*. 
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Proof. Let a = 14 + r" 1 (18 + 25-sJ, and note that since > 22 + 25 + ^_ , r > § , 
we have a < 0. If no is sufficiently large, we may apply proposition 5.2 to obtain, 

n 

|| w n+l 1 1 14 — / II <->iWj ||i4 
i=l 
n 

i=i 

oo 

i=i 

In order to obtain the largest value for s and smallest lower bound for s* which 
satisfy the conditions of propositions 5.1 - 5.4, we choose r = 1.6 so that > 100 
and s < — 96 . We now establish two corollaries which will complete the proof of 
theorem 1.3. 

Corollary 5.1. w n -> w in H s *~ 96 (X 0O ). 
Proof. If s < s* — 96, then by II n 

oo ^ \ V II „, ||fe 



ii w i - w j iir < E ii Uk 

k=j 

< r \^ ll T- 1 (s-s f +18+25) I, , mi 

fc=j 

Hence, {uv} is Cauchy in H s (X oc ) for all s < — 96 since 18 + 25 < 96. □ 
Corollary 5.2. $(w n ) in H 8 *' 96 ^). 
Proof. If s < s* — 96, then by III n 



I 71 

"I lis 



|| *K) L°° < || fn 

< ^r (s - s * +i8+M) iiMii-i 5 ^o. □ 

Since s* > 100, it follows that w n -> u> in C^X^). Therefore $(u> n ) -> 
showing that w is a solution of (5.1). Furthermore, if / is as in theorem 1.3, then we 
have w G C'~ 98 , / > 100. This completes the proof of theorem 1.3. 
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